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SUMMARY 

Ihe  distribution  of  the  point  of  Impact  of  bombs,  missiles, 
and  other  projectiles  frequently  is  approximately  Gaussian 
(normal)  In  character.  If  the  dispersion  is  equal  along  any  two 
perpendicular  axes,  then  the  Impact  probability  law  is  termed 
circular  Gaussian;  If  the  dispersion  Is  not  equal  along  perpen¬ 
dicular  axes,  then  the  probability  law  Is  called  elliptical 
Gaussian.  We  consider  first  the  evaluation  of  the  Integral  of 
an  elliptical  Gaussian  over  a  circle  of  radius  R,  l.e.,  the 
probability  of  a  missile  landing  within  a  circle  of  radius  R 
If  aimed  at  the  center  and  subjected  to  an  elliptical  Gaussian 
Impact  probability  distribution. 

An  equivalent  problem  Is  shown  to  be  the  evaluation  of 
the  Integral  of  a  circular  Gaussian  over  an  ellipse,  Ihls 
Integral,  termed  the  elliptic  coverage  function.  Is  expressed 
as  the  difference  of  two  circular  coverage  functions;  these  have 
been  tabulated  In  RM-330.  A  short  table  of  the  elliptic  coverage 
function  Is  obtained  and  presented  In  a  convenient  graphical 
form.  Finally,  some  characteristics  of  the  circular  probable 
error  and  the  quartlles  of  an  elliptic  Gaussian  are  presented 
In  graphical  form. 
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PREFACE 

Notes  concerning  the  Integral  of  an  elliptic  Oausslan 
distribution  over  a  circle,  taken  by  the  author  from  lectures 
of  H.  H.  Oermond  at  the  University  of  Florida  in  I9^7i  form 
the  basis  of  this  research  memorandum.  The  present  work  extends 
the  results  given  In  RM-330,  The  Circular  Coverage  Function, 
by  H.  H.  Oermond. 

Although  this  work  was  essentially  completed  In  the  early 
1950's,  there  was  Insufficient  Interest  In  the  topic  at  the 
time  to  Justify  the  effort  required  to  put  It  In  shape  for 
external  distribution.  The  past  two  years,  however,  have  seen 
a  general  revival  of  interest  In  Oausslan  Impact  distributions 
of  various  types.  Possibly  this  reflects  a  renewed  attention 
to  methods  for  estimating  the  kill  probabilities  of  hardened 
point  targets,  now  that  city  populations  may  no  longer  be  the 
target  of  Interest.  Since  there  have  been  numerous  requests 
for  any  RAND  work  on  the  topic  of  this  research  memorandum, 
publication  now  seemed  desirable. 


RM-2765 

vll 


co?]TEtrrs 

SUMMARY . Ill 

PREFACE  .  V 

Section 

1.  INTRODUCTION .  1 

2.  THE  CIRCULAR  COVEFJiaE  FUNCTION .  ? 

3.  THE  ELLIPTIC  COVERAGE  FUNCTION .  3 

4.  PROBABILITY  CIRCLES — CEP  AND  CIRCULAR  QUANTITIES.  .  .  13 

Appendix:  THE  CORRELATED  BIVARIATE  DISTRIBUTION . 1? 

REFERENCES . 19 


Rr^2765 

1 


SOME  CHARACTERISTICS  OF  THE  ELLIPTIC  QAUSSIAN  DISTRIBUTION 
1.  INTRODUCTION 

The  distribution  of  the  points  of  Impact  of  bombs,  missiles, 
and  other  projectiles  frequently  is  approximately  Oaussian 
("normal")  in  character.  In  general,  the  dispersion  along  any 
two  mutually  perpendicular  axes  will  not  be  equal  and  the 
correlation  coefficient  will  not  be  zero.  Given  a  map  showing 
the  points  of  Impact,  it  is  convenient  to  use  the  mean  point 
of  Impact  as  the  origin  and  to  choose  the  orientation  of  the 
axes  to  minimize  the  absolute  value  of  the  correlation 
coefficient.  It  is  easily  shown  (see  Appendix)  that  this  choice 
of  orientation  in  the  case  of  a  Gaussian  distribution  is  equlvar- 
lent  to  choosing  axes  oriented  in  the  directions  of  the  maximum 
and  minimum  dispersion,  respectively,  and  that  the  corresponding 
correlation  coefficient  is  zero.  (See  Refs.  1,  2.) 

Without  loss  of  generality,  then,  we  can  thus  restrict 
attention  to  an  uncorrelated  Oaussian  distribution  in  the  xy— 
plane  with  components  of  standard  deviation  a  2 
origin  at  the  mean  Impact  point.  We  designate  such  a  distribu¬ 
tion  as  an  elliptical  Oaussian.  (Equi-probablllty  curves  are 
ellipses  with  major  and  minor  axes  aligned  with  the  x  and  y 
axes,  respectively.)  We  consider  here  the  problem  of  evaluating 
the  Integral  of  such  a  probability  function  over  a  circle  of 
radius  R  (the  probability  of  a  missile  landing  within  a  circle 
of  radius  R  if  aimed  at  the  center  when  the  missile  is  subject 
to  a  Oaussian  impact-probability  law).  An  equivalent  problem  is 
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the  Integration  of  a  symmetric  Qausslan  distribution  over  an 

ellipse.  We  consider  also  the  Inverse  problem  of  determining 

for  a  given  elliptic  Oausslan  distribution  the  radius  of  that 

circle  which  Includes  looq  per  cent  of  the  Impacts  (that  circle 

over  which  the  Integral  of  the  Gaussian  is  q).  Of  particular 

Interest  are  the  0.5  circle,  having  radius  ^  equal  to  the 

"circular  probable  error"  (CEP),  and  the  upper  and  lower  circular 

quartlles,  R_  and  R_  We  consider  the  relation  of  these 

radii  to  the  standard  deviations  o  and  o  and  also  to  the 

A  y 

arithmetic  and  quadratic  means  of  the  o's. 

2.  THE  CIRCULAR  COVERAGE  FUNCTION 

We  Include  here  the  definition  and  some  characteristics  of 
the  circular  coverage  function  that  will  be  needed  later.  In 
Ref.  1,  Germond  designates  the  Integral  of  a  circular  (symmetric) 
Gaussian  distribution  with  unit  standard  deviation  over  a  circle 
of  radius  R  and  at  a  distance  r  from  the  origin  as  p(R,r),  the 
circular  coverage  function.  Thus  p(Ps/a,  r/a)  is  the  probability 
that  a  missile  will  hit  within  a  circle  of  radius  R  If  aimed 
at  a  point  set  off  a  distance  r  from  the  center*  of  the  circle, 
where  the  missile  Is  subject  to  the  Gaussian  Impact-probability 
law  with  zero  correlation  and  symmetric  standard  deviations 
o  -  ■  Cy.  The  function  p(R,r)  Is  expressible  In  the  fol¬ 

lowing  equivalent  f orTma ; 
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(2.1)  p(R,r)  -  exp(-rV2)  J  n  exp(-T)V2)  lo(i’n) 

0 

(2.2)  -  1  -  exp(-rV2)  J  t)  exp{-^^/2)  lQ{rr\)  d^ 

R 

(2.3)  «  exp[-(R2  +  r^)/2)  S  dl/r)*^ 

nol  ^ 

(2.4)  -  1  -  expHR^  +  r^)/2)  Z  (r/R)*^ 

n»0 

where  lyj(z)  is  the  modified  Bessel  function  of  the  first  kind, 
of  order  n.  Tables  of  p(R,r)  have  been  prepared  In  Ref.  1. 

3.  THE  ELLIPTIC  COVERAGE  FUNCTION 

We  designate  the  Integral  of  a  symmetric  oaussian  distri¬ 
bution  of  unit  standard  deviation  and  zero  correlation  co¬ 
efficient  over  an  ellipse  with  major  axis  a  and  minor  axis  p 
as  q(a,p),  the  "elliptic  coverage  function."  Thus  q(a/a,  p/a) 
Is  the  probability  that  a  missile  will  Impact  within  the 
ellipse  (a,p)  If  aimed  at  the  center  when  the  missile  Is 
subject  to  a  aausslan  impact-probability  law  with  zero  corre¬ 
lation  and  symmetric  standard  deviations  a.  We  show  that  It 
Is  expressible  as  the  difference  of  two  circular  coverage 
functions,  and  Is  thus  also  a  tabulated  function.  Further, 
we  show  that  the  Integration  of  an  elliptic  Oaussian  over  a 
circle  Is  an  equivalent  problem.  Finally,  using  the  tabulated 
values  of  the  circular  coverage  function,  we  present  a  simpli¬ 
fied  graphical  representation  of  the  function  q(lV'a  ,  B/a  ). 

y  ^ 
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We  first  prove  three  lemmas  that  will  later  be  needed. 
Lemma  1 .  For  p,  real^  ||i|  <  1,  we  have 

1  —  (JL^  00 

(3.1)  - n - =  1+2  S  (-1)"^  M-”  cos  (nq)). 

1  +  +  2ii  cos  <p  n=rl 


Proof.  Consider  the  congjlex  variable  z  ■  Xe^*^.  For 
Tv  <  1,  the  function  (1  —  z)/(l  +  z)  has  the  power  series 
expansion 


1  —  z 
1  +  z 


1  +  2  E  (-1)*^  . 

n=l 


Taking  real  values  of  both  sides,  we  obtain 

1  -  Ts^  00 

- n -  =  1+2  S  (-1)’^  Tv”  cos  (ncp). 

1  +  Tv*^  +  2X  cos  cp  n«l 


Hius  the  lemma  holds  for  p.  positive.  For  p.  negative,  the 
proof  is  similar,  involving  now  (1  +  z)/(l  -  z)  in  place  of 
(1  -  z)/(l  +  z). 

Z  Cos  CD 

Lemma  2.  Hie  function  e  ^  for  all  z  and  <p  may  be 

expressed  by  the  expansion 

(3.2)  e^  ^  =  Iq(z)  +  2  Z  i^(z)  cos  (nip), 

n=l 

where  is  the  modified  Bessel  function  of  the  first  kind,  of 
order  n. 

Proof.  Ihls  well-known  expansion  in  Bessel  functions  is 
most  easily  proved  in  terms  of  the  Bessel  coefficients.  It 
may  be  verified  from  the  definitions  of  Bessel  functions  that. 
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for  all  z  and  nonzero  values  of  t, 

”*■  ^  ^  * 

Letting  t  e  we  see  that 

W  j 

exp(z  cos  (p)  *1  S  I^(z)e^ 
n»-«o 

Equating  real  parts  and  noting  that  I_^(2)  « 

cos  <P  „  Iq(z)  +22  * 

nal 

Lemma  3 .  If  p(R,r)  Is  the  circular  coverage  function  of 
Sec.  2j  then 

(3.3)  p(R,r)  -  p(r,R)  -  1  -  exp[~  (R^  +  r^)/2] 

,  rio(rR)  +  2  2  In(Rr)l- 

i_  n«»l  J 

Proof.  Using  the  form  (2.3)  for  p(r,R)  and  (2.4)  for 
p(R,r)>  we  obtain 

00 

p(R,r)  -  p(r,R)  -  1  -  exp[-.(R2  +  r^)/2]  2  (|)  I  (Rr) 

n-0  ” 

-  exp[-  (R^  +  r^)/2)  2  (J)*^ 

ns»l 

whence  (3.3)  follows. 

Theorem  1.  The  elliptic  coverage  fxmctlon  q(A,B),  the 
integral  of  a  symmetric  Gaussian  distribution  of  unit  standard 
deviation  over  an  ellipse  with  major  axis  A  and  minor  axis  B, 

Is  given  by  the  difference  of  two  circular  coverage  functions 
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(3.4)  q(A,B)  .  p(A^,i_^)  _ 


Proof.  By  definition. 


S  +  y^)/2a^]dx  dy, 

Stto 


where  the  integrations  are  over  the  area  bounded  by  the  curve 
2  2  2  2 

C.*  X  /a  +  y  /p  n  1.  Letting  A  «  a/a  and  B  »  p/a, ^  ■  x/a  and 
T)  »  y/o,  we  obtain 

(3.5)  q(A,B)  expl-  +  n^)/2]  d^  dii 

C'^ 

where  C ' :  $ «  1 . 

Making  a  transformation  to  polar  coordinates,  $  -  p  cos  B, 
T]  =*  p  sin  e,  and  letting 


/ cos^  6 

^ir~ 


+ 


sin^ 


from  (3.5)  we  obtain 

(3.6)  q(A,B)  f  P  exp(-  p^/2)  dp  d© 


2Tr  L 

I 

0  0 


2Tr 


Making  the  further  transformation  tan  6  -  ^  tan  <p,  we  have 

2  2  2  2  2 
L  =  A  cos  <p  +  B  sin  <p, 

de  -  AB  d(p/L^, 


and  thus 
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(3.7)  q(A,B)«l-4/ 


IT  p 

e:3ro(-  L  /2)AB  dq 


»2  .  „2i 


.2  „2. 


coa^  9  +  B^  Bln^  9  ••  t  +  lA — cos  2q), 


BO  (3.7)  becomes 

(3.8)  q(A,B)  -  1  -  expf-  ( 


•  I  ~  /  exp  [(-  A-.  cos  29]  — .  . 

I  0  ^  - 


Further,  letting  p.  «  (A  —  B)/(A  +  B)*  we  have 


L  A  +  B  +  (A  ~  B  )cos  29  1  +  +  2m.  cos  29 

Utilizing  Lemma  1,  we  can  write 

(3.9)  ^  o  1  +  2  2  (-1)^  p.’^  COB  2n9 

L  n=.l 

“1  +  2  2  (—1)'^  • 
n=l  a  + 

From  Lemma  2  we  obtain  the  expansion 

(3.10)  exp[-  cos  29]  -  Iq(A...~..B^) 


n  T  /A^  -  eP 


2  2  (-1)“ 
neil 

Substituting  (3.9)  and  (3.10)  In  (3.8),  and  using  the 


)cos  2n9. 


r  IT 

/  cos  2m9  cos  2n9  d9 


fact  that 
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for  m  o  n  and  zero  otherwise,  we  obtain 


(3.11)  q(A,B)  =  1  -  e 


If  we  let 


A  +  B  n  A  —  B  j, 

n  “  n,  A  ^  t 


It  follows  that 


2  2  2  2  2  2 
il+^„£4X-  and  Hr, 


SO  that  (3.11)  becomes 

(3.12)  q(A,B)  -  1  ~  exp(-  i^(Rr)  +  2^^2^(5)”  X^irR)  • 

By  Lemma  3  the  right-hand  side  of  (3.12)  Is  p(R,r)  —  p(r,R), 

BO 

q(A,B)  »  p(R,r)  -  p(r,R)  «  p(^A_+_S, 

Ya  -  B  A  +  B^ 

-  2  '  5  7* 

Theorem  2.  The  integral  of  an  elliptical  Gaussian  distrl— 

butlon  with  components  of  standard  deviation  a„  >  a„  and 
"  X  y 

origin  zero  over  a  circle  of  radius  R  is  given  by  the 


elliptical  coverage  function 


■<t'  4)- 


Proof.  The  Integral  J  of  the  elliptical  Gaussian  dis¬ 
tribution  over  a  circle  of  radius  R  is  given  by 
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^  -  h  Vy //  j-  dy, 

0^0 

where  Ci  x  +  y  «  r  ♦  under  the  traneformatlon  x/ex  « 
y/cry  »  t),  we  obtain 

J  -  exp[l(^^  +  Ti^)/2]d§  dT) 

C 

where  O'  i  ox^^^/r^  +  oy^T)^/R^  •»  1.  Thue,  by  definition  of  the 
elliptic  coverage  function, 

'  ■  $■ 

Using  Hieorem  1  and  the  tabulated  values  of  the  circular 
coverage  function  of  Ref.  1,  we  present  in  Table  1  a  short 
table  of  values  of  the  elliptic  coverage  function  q!^~# 

For  convenience  we  have  used  as  entries  the  ratios 
p  «  where  Is  the  maximum  standard 

deviation.  Thus 


$  ■  45;. 


1  “  R  A  +  p'\ 


2p 


r-r  2p 


For  p  «  0,  l.e.,  o  =  0,  the  distribution  q  becomes  the 

O' 


linear  Qausslan,  l.e., 

11m  /-H,  jl)  .js  r 

a„-40^VV 


^x  p 

exp(-x  /2)  dx. 
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The  data  of  Table  1  are  presented  In  a  very  convenient 
graphical  form  In  Fig.  1. 


6  66 
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R/tr„ 

Fig.  I  —  Elliptic  coverage  function 
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4.  PROBABILITY  CIRCLES — CEP  AND  CIRCULAR  QUANTITIES 

In  Sec.  3  we  have  determined  the  probability  q  in  terms 

of  the  circle  radius  R  and  the  standard  deviations  a_  and  a 

A  y 

It  Is  often  desirable  to  determine,  for  a  fixed  value  of  q, 
the  radius  of  that  circle  for  which  the  elliptic  coverage 
function  (probability  of  landing  within  R^)  is  q.  The  more 
commonly  used  of  these  are  the  0.5  circle,  whose  radius  R^  ^ 
is  called  the  CEP,  and  the  0.25  and  0.75  circles  whose  radii 
we  will  call  th*e  lower  and  upper  circular  quartlles,  respec¬ 
tively. 

Figure  1,  which  shows  q  as  a  function  of  p  and 
can  also  be  considered  to  show  as  a  function  of  q  and  p. 

Let  h(q,p)  Then 

Rq  -  # 


where  h  is  graphed  in  Pig.  1.  It  is  often  desirable  to  relate 


Rq  to  the  arithmetic  or  quadratic  means  of  the  o's,  l.e.. 


o  +0 

— y 


2  ,  2 
2  °y 

'  2" 


In  Pig.  2  the  ratios  of  Rq/o  for  q  »  0.25,  0.50,  0.75  are 
plotted  against  the  ratio  p,  while  in  Fig.  3  the  same  is  done 
for  the  ratio  R  /q. 

q 

A  particular  advantage  of  the  ratio  Rq/Q  is  that  the 
quadratic  mean  Q  (see  Appendix)  is  Independent  of  correlation. 
Of  particular  interest  in  Pigs.  2  and  3  are  the  flatness  of  the 
curves  and  R^  l.e.,  for  q  of  the  order  of  0.5, 
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of  the  order  of  0-75*  VQ  independent  of  p  for  all  p 
we  have  as  approximations 

Rq  ^  ®  CHP  1«17 

p 

0.75  o  upper  circular  quartlle  »  1.65  Q. 


IhUB 
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Fig.  3 


o  ..  Minimum  <r 

Ratio,— — : -  =  p 

Moximum  <r 


Ratio 


as  0  function  of  p 
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Appendix 

THE  CORRELATED  BIVARIATE  DISTRIBUTION 

Consider  a  bivariate  Oausslan  distribution  In  the 

xy— plane,  with  center  at  the  origin,  with  components  of  standard 

deviation  equal  to  o  and  o  ,  and  with  a  coefficient  of 

X  y 

correlation  equal  to  r.  The  equation  of  such  a  distribution 
can  be  written 

dp - ^ - -  expf-  iLj/g  _  r®)  . 

2TO^OyNl-r  L  Vy  “y  -I 


Fig.  4 


The  substitution 


X  -  u 

cos 

CD 

1 

< 

sin 

e 

y  »  u 

sin 

e  +  V 

cos 

Q 

where 
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tan  20  « 

rotates  the  axes  and  converts  this  Into  an  uncorrelated, 
bivariate  Gaussian  distribution: 


Vv  ,u2  ^  v^ 

U  V  [ 


The  components,  and  o^,  of  standard  deviation  along  the 

new  axes  are  given  by 


and 


regardless  of  the  value  of  0  (and  hence  r.-igardless  of  the 
value  of  r ) . 
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